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ABSTRACT 


Multi-target  tracking  ( MTT )  has  many  applications,  and  has 
therefore  been  the  subject  of  considerable  investigation.  One  key 
aspect  of  this  problem,  especially  in  a  dense  target  environment, 
is  the  "scan-to-scan  correlation"  or  "tracking  data  association" 
problem  of  assigning  measurements  to  tracks.  Various  different 
approaches  have  been  proposed  to  solve  this  problem.  In  this 
report,  three  algorithms  for  MTT  data  association  are  presented  for 
comparison.  First  the  nearest-neighbor  standard  filter  (NNSF) 
algorithm  is  presented.  Then  two  of  the  more  promissing  extensions 
are  presented:  the  multiple  hypothesis  test  method  ( MHT ) ,  and  joint 
probabilistic  data  association  (JPDA)  method.  Whenever  possible  the 
same  notation  is  used  in  presenting  all  three  methods  for  ease  of 
comparison.  This  report  is  intended  to  serve  as  a  prelude  to  a 
comparative  investigation  of  these  three  competing  data  association 
methods . 


RESUME 


Les  techniques  de  poursuites  de  cibles  multiples  ont  plusieurs 
applications  et  ont  done  fait  l'objet  de  nombreuses  recherches. 
L'un  des  aspect  important  du  probleme  est  1 r association  des  donnees 
radar  avec  les  pistes.  Plusieurs  approches  ont  deja  ete  proposees 
pour  resoudre  ce  probleme.  Ce  rapport  compare  trois  algorithmer. 
d '  association  des  donnees.  La  methode  standard  du  voisin-immediat 
(NNSF)  est  d'abord  presentee.  Deux  variations  prometteuse  sont 
ensuite  analysees:  une  methode  de  test  a  hypotheses  multiples 
(MHT)  et  une  methode  d' association  de  donnees  probalistique  (JPDA). 
Ce  rapport  presente  une  analyse  preliminaire  des  trois  methode s  et 
sera  suivi  d'etudes  plus  approf ondies . 
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EXECUTIVE  SUMMARY 


In  this  report  three  quite  different  algorithms  for  multiple 
target  tracking  in  a  cluttered  environment  are  presented:  Nearest 
Neighbour  Standard  Filter  ( NNSF ) ,  Multiple  Hypothesis  Test  (MHT) 
and  Joint  Probabilistic  Data  Association  (JPDA).  The  basic  tracking 
equations  and  data  association  algorithms  are  presented,  using  a 
consistent  notation,  so  that  the  different  methods  can  be  directly 
compared.  The  apparent  advantages  and  dissadvantages  of  each  are 
briefly  explored,  and  suggestions  made  for  further  analysis. 

Superficially  the  tradeoff  among  the  three  different 
approaches  is  simply  one  of  accuracy  vs.  computational  efficiency, 
with  the  NNSF  method  being  the  least  accurate  but  most  efficient, 
and  the  MHT  method  being  the  most  accurate  and  least  efficient. 
However  this  simplistic  view  overlooks  the  fact  that  the  three 
methods  differ  fundamentally  in  their  approach  to  data  association. 

The  NNSF  produces  a  single  unambiguous  data  association 
solution  at  each  point  in  time,  based  on  the  previous  association 
and  the  current  sensor  information.  This  however  may  not  be  the 
best  choice,  especially  because  it  does  not  make  full  use  of  all 
prior  sensor  data.  Once  an  incorrect  association  is  made,  it  seems 
unlikely  that  the  solution  would  ever  recover. 

The  MHT  approach  will  maintain  several  (perhaps  many) 
different  possible  data  association  solutions,  and  uses  the  history 
of  sensor  data  to  eliminate  highly  unlikely  choices,  eventually 
leaving  only  one  best  choice  (hopefully).  This  should  yield  the 
best  solution,  however  it  does  generally  have  periods  of 
uncertainty.  Thus  if  one  could  not  wait  for  the  solution,  then  it 
would  be  necessary  to  take  special  measures. 

The  JPDA  approach  does  not  explicitly  make  a  data  association 
decision.  Instead  this  method  applies  all  ambiguous  measurements 
to  each  tracking  filter  with  which  it  could  reasonably  be 
associated  (with  an  appropriate  weighting  factor).  Thus  any  given 
tracking  filter  will  likely  be  assigned  several  "incorrect" 
measurements  (hopefully  with  a  low  weighting)  but  will  almost 
certainly  also  be  assigned  the  "correct"  measurement  (hopefully 
with  a  high  weighting).  It  is  thus  hoped  that  the  incorrect 
measurements  (hopefully  random)  will  have  a  small  cumulative 
effect . 

Further  investigation  would  be  necessary  to  determine  which, 
if  any,  of  these  data  association  methods  is  most  appropriate  for 
a  given  multi-target  tracking  problem. 
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1.0  INTRODUCTION  TO  MULTIPLE  TARGET  TRACKING 


1.1  The  Problem 


The  multiple  target  tracking  problem  is  encountered  n  many 
situations,  whenever  sensor  data  is  available  from  one  or  more  wide 
angle  of  view  sensors  (such  as  radar,  electronic  support  measures 
(ESM),  infra-red  (IR),  video,  etc,),  providing  information  on  the 
position  (such  as  range,  azimuth,  elevation  etc.)  of  "targets'' 
(aircraft,  vessels,  missiles,  ground  vehicles  etc.).  TS.e  sensors 
may  also  provide  information  useful  for  target  identification 
purposes  (size,  velocity,  luminosity,  transmitter  frequency,  etc.). 
Although  this  information  may  be  useful  to  a  target  tracking 
algorithm,  the  identification  aspects  of  target  tracking  are  beyond 
the  scope  of  this  report. 


These  target  tracking  sensors  typically  scan  their  field 
of  view  at  regular  intervals.  There  may  be  many  contacts  from  such 
a  scan,  from  the  different  objects  in  the  sensors  field  of  view, 
as  wTell  as  some  which  may  be  just  noise.  The  sensor  measurements 
from  all  contacts  within  a  scan  are  generally  processed  together, 
since  they  can  safely  be  assumed  to  correspond  to  different  objects 
(something  which  cannot  be  said  of  returns  from  different  scans  or 
from  different  sensors). 


The  basic  tracking  problem  is  to  estimate  the  position  and 
velocity  of  the  target(s),  using  the  available  sensor  data  (from 


a  sequence  of  scans), 
predict  the  targets' 
is  a  problem  because 
ambiguities  and  the 
prediction  is  further 


It  is  also  normally  desireable  to  be  able  to 
location  some  time  in  the  near  future.  This 
the  sensor  data  normally  has  errors  and/or 
targets  generally  move  between  scans.  The 
complicated  by  unknown  target  manoeuvres.  If 
there  is  only  a  single  target,  however,  then  Kalman  filtering 
techniques  provide  an  "optimal”  solution.  Although  the 
nonlinearitites  of  this  problem  introduce  some  complications  in  the 
use  of  Kalman  filtering,  this  problem  has  been  extensively  studied 
and  suitable  solutions  exist. 


The  multi-target  tracking  problem  however,  is  not  simply 
a  tracking  problem  when  there  are  more  than  one  target.  As  shall 
be  seen  below  however,  the  possibility  of  more  than  one  target 
significantly  complicates  the  situation,  because  it  introduces  the 
problem  of  associating  measurements  with  targets,  when  the  distance 
between  targets  is  comparable  to  the  measurement  error  or  to  the 
error  in  target  position  prediction  (from  one  scan  period  to  the 
next),  then  this  association  problem  can  be  very  difficult.  For 
example  fighter  aircraft  flying  in  tight  formation  being  tracked 
by  a  scanning  radar. 
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1.2  Background 


There  currently  is  no  satisfactory  unified  solution  method 
for  the  Multi-target  tracking  problem.  The  single  target  tracking 
problem,  however,  has  been  studied  in  depth,  and  optimal  estimation 
theory  (Kalman  filter  theory  etc.),  does  provide  a  satisfactory 
solution  (see  for  example  [1)).  With  this  method,  one  can  form  the 
statistically  optimal  estimate  of  the  target  state  vector  (position 
and  velocity)  by  recursively  processing  sensor  measurements  (radar,, 
sonar,  ESM  etc.)  taken  from  the  target  being  tracked. 

The  multi-target  tracking  problem  can  be  approached  in  the 
same  way,  as  a  number  of  stochastic  estimation  problems  (one  for 
each  target).  Unfortunately  there  is  the  added  difficulty  of 
deciding  which  measurements  correspond  to  which  targets  (and  which 
measurements  don't  correspond  to  any  targets).  This  problem  of 
associating  measurements  with  targets  is  known  as  the  "data 
association  problem".  Since  this  is  not  a  problem  of  estimating  the 
value  of  a  random  variable,  it  does  not  fall  neatly  into  the  realm 
of  optimal  estimation.  However,  it  is  amenable  to  stochastic 
treatment,  and  several  approaches  have  been  suggested.  This  report 
describes  in  brief  three  of  the  more  promising  approaches. 

All  three  approaches  are  based  on  the  intuitive  premise 
that  the  probability  of  a  particular  measurement  corresponding  to 
a  specific  target  is  a  function  of  the  proximity  of  that 
measurement  to  the  expected  target  location. 

The  simplest  of  the  approaches  described  here  is  the 
"Nearest  Neighbor"  method,  which  at  each  time  epoch  (scan)  makes 
a  complete  assignment  decision.  This  complete  assignment  can  be 
defined  as  a  specific  decision  for  each  measurement: 

1-  that  it  came  from  one  of  the  targets  being  tracked, 

2-  that  it  is  noise,  or 

3-  that  it  is  from  a  new  target. 

Furthermore  this  allocation  of  measurements  is  made  in  such  a  way 
that  no  target  is  assigned  more  than  one  measurement  from  any  given 
sensor  scan.  The  nearest  neighbor  approach  makes  this  complete 
assignment  decision  at  each  point  in  time  by  minimizing  a  global 
distance  function  which  represents  the  closeness  of  each 
measurement  to  its  assigned  target.  A  problem  with  this  approach 
is  that  if  an  incorrect  assignment  is  made  for  a  single 
measurement,  then  the  tracking  filter  that  processes  it  will  make 
a  poor  prediction  for  the  position  of  it's  target  for  the  next  time 
epoch,  likely  leading  to  subsequent  incorrect  assignments,  and  thus 
the  process  could  easily  break  down.  The  two  other  methods 
described  here  are  different  approaches  to  solving  this  problem. 
The  nearest  neighbour  method  is  described  in  more  detail  in  chapter 
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2  below. 


The  second  approach  to  be  described  is  the  Multiple 
Hypothesis  approach.  As  the  name  implies,  this  approach  attempts 
to  solve  the  problem  (of  not  being  able  to  correctly  decide  at  each 
epoch  which  measurement  corresponds  to  which  target)  by  keeping  all 
reasonably  likely  "complete  assignments"  (as  defined  above)  as 
working  hypothesis  (with  which  to  predict  the  next  epoch's  target 
positions).  The  global  distance  measure  (of  proximity  of 
measurements  to  predicted  targets)  can  then  be  used  to  assign  a 
"probability"  to  each  hypothesis  at  each  epoch.  These  probabilities 
can  then  be  combined  over  several  epochs,  with  the  expectation  that 
incorrect  hypothesis  will  lead  to  highly  unlikely  cumulative 
probabilities  (ie.  dead  ends).  In  this  way  these  hypothesis  can  be 
dropped  and  the  correct  hypothesis  should  eventually  stand  out 
alone.  The  problem  with  this  method  is  that,  unless  the  number  of 
hypothesis  carried  forward  from  one  epoch  to  the  next  is  kept  very 
low,  then  the  "hypothesis  tree"  will  grow  extremely  quickly,  and 
a  computationally  unmanageable  situation  will  arrise.  This  method 
is  described  in  more  detail  in  chapter  3  below. 

The  third  method  is  the  "Joint  Probabilistic  Data 
Association"  or  JPDA  method,  which  takes  a  very  different  approach: 
assigning  all  measurements  in  the  vicinity  of  a  given  target  to 
that  target,  but  weighting  the  me a sue  emeu Is  according  to  their 
proximity.  Thus  it  in  highly  unlikely  that  the  correct  measurement 
will  be  ignored,  but  on  the  other  hand  it  is  an  almost  certainty 
that  incorrect  measurements  will  also  have  some  unwanted  influence. 
This  JPDA  method  is  described  further  in  chapter  4  below. 


1.3  State  Space  Description 


For  an.  object  being  tracked,  the  discretized  equations  of 
motion  may  be  modeled  by 


xk,i  =  Mk  +  Gkuk 


(1) 


where  xk  is  an  nxl  state  vector  (describing  the  position  and 
velocity  etc.)  of  the  tracked  objects  at  the  kth  sample  time, 
is  the  transition  matrix,  and  uk  is  a  mxl  state  excitation  vector 
to  account  for  both  maneuvres  and  modeling  errors  and  is  generally 
assumed  to  be  white  and  Gaussian  with  zero  mean  and  covariance  Qk . 
In  a  track-while-scan  system,  the  kth  sample  will  occur 

approximately  at  time  kT,  where  T  is  the  scan  interval  of  the 
sensor . 


The  observation  equation  representing  valid  sensor 


3 


measurements  of  the  objects  being  tracked  has  the  form 


y*  -  Hkx.<  +  vk 


(2) 


where  yk  is  the  mxl  sensor  measurement  vector  (also  referred  to  as 
the  sensor  "reports")  and  vk  is  white  Gaussian  observation  noise 
with  zero  mean  and  covariance  Rk .  The  observation  equation  for 
extraneous  sensor  reports  resulting  from  thermal  false  alarms, 
clutter,  and  other  targets  is  assumed  to  satisfy 


Yk  -  HkXe^^j  +  wk 


(3) 


where  wk  is  assumed  white  and  uniformly  distributed  over  some 
volume  C  of  the  measurement  space  centered  about  the  predicted 
measurement 


Yek  |  k-1  *  Hkxek|k-1 


(  3a) 


The  number  of  such  extraneous  reports  in  any  volume  C  obeys 
a  Tcisson  distribution  with  mean  BC  ,  where  B  is  the  unnormalized 
extraneous  reports  density. 

The  tracking  filter  provides  a  state  estimate  xek|k  and 
one-scan-predicted  state  xek+1,,(  given  all  measurement  data  up  to 
the  time  k.  The  basic  filter  equations  are  described  in  the 
following  section. 


14  Basic  Tracking  Filter  Equations 


The  tracking  filter  state  estimate  xek|k  and  covariance 
Pk|k  update  equations  are  calculated  (as  functions  of  the 
measurements  yk  and  the  previous  estimates  xek|k_x  and  Pk|k_1)  as 
follows: 


xeMk 

“  xek|k-l  +  ^k^Yk  -  ^kXek|k-l} 

(4) 

it 

K 

^k|k-l^k  l  ^k^klk-l^k  +  Rk  J 

(5) 

Pk|k 

■=  [  I  -  KkHkJPk|k_1 

(6) 
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This  requires  knowledge  of  the  measurement  noise  covariance  Rk  and 
the  measurement  matrix  Hk  relating  the  measurements  to  the  state 
vector . 


These  estimates  can  then  be  propagated  one  scan  into  the 
future  (prediction)  by  using  the  state  transition  matrix  as 
follows : 


xek<i|k  K  *kxek|k  n> 

^k  +  llk  “  *kPk|k^k  +  ^kQkGk 


1.5  Notation  and  Definitions 


To  facilitate  comparison,  the  same  notation  is  used 
wherever  possible  to  describe  the  different  approaches  to  data 
association.  This  notation,  and  related  definitions  are  collected 
here  for  ease  of  reference: 


k  =  time  or  scan  index. 

N [ x , P ]  =  probability  density  function  for  a  zero  mean  Normally 

distributed  (Gaussian)  random  vector  x,  with  covariance 
matrix  P 

-  exp(-0.5  xT  P"1  x)//{(2rc)n  |  P  |  }  (9) 

(where  n  =  dimension  of  x). 

Nc  <=  no.  of  measurements  associated  with  confirmed  tracks. 

Nf  »  number  of  measurements  associated  with  false  targets. 

Nn  -  number  of  measurements  associated  with  new  targets. 

Nt  *=  number  of  previously  known  (at  time  k-1)  targets 
(confirmed  tracks)  within  the  cluster. 

Nm;  number  of  measurements. 


Pp  -  probability  of  detection. 
Bft  =  density  of  false  targets. 
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Skt  «*  density  of  previously  known  targets  that  have  been 
detected . 

c  «  normalization  constant  given  by  summation  of  all  the 

hypothesis  probabilities  within  the  cluster  for  one  scan, 

Pr(k)  «  probability  of  a  hypothesis  at  time  k  given  probability  of 
the  prior  hypothesis. 

Zb„1L  ■=  Summation  over  b  from  1  to  L. 

lib.!1,  -  Product  over  b  from  1  to  t . 

Lk  «  total  number  of  hypothesis  for  track  1  at  time  k  in  the 
absence  of  pruning: 

Lk  =  (  1  +Nk  )  Lk_j 

for  N'k  the  number  of  sensor  reports  falling  within  the  gate 
for  track  i  at  scan  k  (pruning  may  have  reduced  Lk_1 
between  scans ) . 

Palk  |  k_i  =  covariance  of  the  estimation  error  given  observations 
through  scan  k-1  and  given  track  hypothesis  "  a  1 "  . 

Ka 1 k  j  v  - 1  =  expected  value  of  the  estimation  error  given  observations 
through  scan  k-1  and  given  track  hypothesis  "al". 

Pak|k_1  and  Kak|k-i  are  the  an<3  Kal  analogues  corresponding  to 

the  hypothesis  "a". 

S_1k  =  inverse  of  Sk  which  is  the  covariance  of  the 

innovation  yk  -  ye^.!  given  hypothesis  "al". 

xeMk-i  “estimated  state  at  scan  k  given  hypothesis  up  to  time  k-1. 

pk|k  =  covariance  of  the  estimation  error  given  observations 
up  to  scan  k. 

Ak  =  optimal  tracking  filter  correction  vector  for  MHT  filter. 

AAk  *=  optimal  tracking  filter  gain  for  JPDA. 

De  -  the  expected  track  length. 

=  the  event  that  an  observation  j  belongs  to  a  target  t. 

Zk  =  the  set  of  measurements  up  to  time  k. 

Pp1  =  the  probability  of  detection  of  target  t. 

N t ( j )  =  number  of  targets  associated  with  the  observation  j 
except  for  false  alarms(  false  targets). 
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Nc(t)  -  number  of  observations  associated  with  the  target  t 

Nc'(t)  "  number  of  observations  not  associated  with  the  ta.get 

Ncr  «*  (Number  of  columns  of  9)*(Number  of  rows  of  2), 
where  2  is  validation  matrix  given  by  eq.(25). 


2.0  A  NEAREST  NEIGHBOR  STANDARD  FILTER  ALGORITHM 


2.1  introduction 


The  nearest  neighbor  approach  to  data  association 
determines  a  unique  pairing  so  that  at  most  one  observation  from 
each  sensor  scan  can  be  paired  with  a  previously  established  trac::. 
Also,  a  given  observation  can  be  used  only  once,  either  to  update 
an  existing  track  or  to  start  a  new  track  (which  may  be  abandoned 
os  noise  if  no  subsequent  observations  match  it).  This  method  is 
based  on  likelihood  theory,  and  the  goal  is  to  minimize  an  overall 
distance  function  that  considers  all  observation-to-track  pairings 
that  satisfy  a  preliminary  gating  test. 


2.2  NNSF  Algorithm 

In  this  section,  a  NNSF  algorithm,  based  on  the  modified 
Munkres  optimal  assignment  method  (see  reference  [2]),  is 
presented . 


Step  1:  Set  k*l  (scan  index)  and  all  control  parameters. 

Step  2:  Initialize  the  sy~tem  parameters;  *,  Q,  H,  R. 

Step  2:  Simulate  target  trajectories. 

Step  4:  Simulate  measurements  (or  receive  the  data  from  radar.) 

Step  S:  Validate  the  measurements  using  gating  test:  Form  gate 
about 

yekik-!  *  H  *  xek|k_! 

and  select  the  Nk  sensor  reports  to  be  used  in  filter 
updating.  A  measurement  is  valid  if 

{yk  -  yeMk-i)T  S'\  (Yk  -  YeMk-i)  ^  9*  (10) 

(this  is  a  'g-o"  elliptical  gate). 

Step  6:  Form  the  assignment  matrix.  The  elements  of  the  matrix 

are  equal  to  the  normalized  distance  function  associated 
with  the  assignment  of  each  of  "Nobs"  observations  to  each 
of  "Ntr"  tracks.  If  the  gating  relationship  is  not 
satisfied,  the  observation-  to-track  pairing  can  be  given 
a  very  large  distance  to  penalize  this  assignment. 


8 


Step  7:  Solve  the  Assignment  Matrix:  Minimize  the  normalized 

distance  function  using  the  modified  Munkres  optimal 
assignment  algorithm  in  section  2.3. 

Step  8:  Correlate  the  observations  to  the  tracks  according  to  the 
solution  of  optimal  assignment  matrix. 

Step  9:  Update  and  predict  the  state  vectors  (tracks)  using  (4) 
to  (9). 

Step  10:  Output  the  predicted  track  positions. 

Step  11:  If  k=kf  (end  of  tracking  mission),  go  to  the  Step  11. 
Otherwise,  set  k-k+1  and  go  to  Step  4. 


2.3  Modified  Munxers  Assignment  Algorithm 


The  following  method  is  based  on  the  Munkres  optimal 
assignment  algorithm  modified  by  Burqeois  and  Lassalle  [12],  This 
method  has  an  advantage  (compared  to  older  methods)  for 
applications  since  the  assignment  matrix  need  not  be  square.  For 
the  convenience  of  presentation  of  the  algorithm,  the  rows  and 
columns  of  the  matrix  may  be  marked  and  referred  to  as  covered.  The 
zeros  may  be  marked  by  being  starred  (*)  or  primed  ('). 


The  Optimal  Assignment  Algorithm: 


Step  1: 

Step  2: 
Step  3: 
Step  4: 

Step  5: 


Initially,  no  lines  are  covered  and  0's  are  s-arred  or 
primed. 

Let  v  *  min. {Number  of  rows,  Number  ot  columns}. 

If  (number  of  rows)  >  (number  of  columns),  go  to  Step  6. 

For  each  row  in  the  matrix,  subtract  the  smallest  element 
of  the  row  from  each  element  in  the  row. 

If  the  (number  of  columns)  >  (number  of  rows),  go  to  Step 
7  . 


Step  6:  For  each  column  in  the  matrix,  subtract  the  smallest 
element  of  the  column  from  each  component  of  the  column  and 
go  to  Step  7 . 

Step  7:  a)  Find  a  zero,  "Z",  of  the  matrix. 

b)  If  there  is  no  starred  zero  in  its  row  or  its  column, 
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star  the  zero  (i.e.  Z*).  Repeat  for  all  zeros  of  the 
matrix.  Go  to  Step  8. 

Step  8:  a)  Cover  every  column  containing  a  starred  zero  Z*. 

b)  If  v  columns  are  covered,  the  locations  of  the  Z*  form 
the  row-column  associations  (i.e,  observation-to-t.rack 
pairs).  The  algorithm  is  now  completed.  Otherwise, 
continue  to  the  next  step. 

Step  9:  a)  Choose  an  uncovered  zero  and  prime  it  (i.e.  Z'). 

b)  If  there  is  no  starred  0  in  the  row  of  Z',  go  to  Step 

10. 

c)  If  there  is  a  starred  zero  Z*  in  the  row  of  Z',  cover 
this  row  and  uncover  the  column  of  Z*. 

d)  Repeat  until  all  zeros  are  covered  and  then  go  to  Step 

11. 

Step  10:  al)  Let  Zo  denote  the  uncovered  2'.  If  there  is  no  Z*  in 
the  column  of  Zo,  go  to  Step  (a6). 
a2)  Let  Z 1  denote  the  Z*  in  the  column  of  Zo. 
a3)  Let  Z2  denote  the  Z'  in  the  row  of  Zl. 
a4)  Continue  the  steps  a2  and  a3  until  a  7,2  which  has  no 
Z*  in  its  columns  has  been  found. 
a5)  Un-star  each  starred  zero  of  the  sequence, 
a  6 )  Star  each  primed  zero  of  the  sequence, 

bl)  Erase  all  primes  from  primed  zeros  and  uncover  every 
line . 

cl)  Go  to  Step  8. 

Step  11:  a)  Find  the  smallest  uncovered  element  in  the  matrix  and 
call  it  "m".  "m"  will  be  positive. 

b)  Add  "m"  to  each  covered  row. 

c)  Subtract  "m"  from  each  uncovered  column. 

d)  Go  to  Step  9  without  altering  stars,  primes,  or 
uncovered  lines. 
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3.0  THE  MULTIPLE  HYPOTHESIS  APPROACH 


3.1  Introduction 


In  the  sequel  of  this  section,  a  track  history  "a"  at  scan 
k  is  defined  by  selecting  a  single  sensor  report  y  from  each  scan 
j<k 


j);  j  =  1 , 2 .  . . k  I  0  <1 j  <N j } 


where  Nj  is  the  number  of  reports  at  scan  j  and  lj  =  0  refers  to  the 
hypothesis  that  none  of  the  sensor  reports  within  the  gate 
originated  from  the  target.  Hence  the  track  history  "a"  is  just  the 
hypothesis  that  the  entire  sequence  of  measurements  within  "a"  is 
correct;  i.e.  each  sensor  report  y(Ij;j)  originated  from  the  target 
when  Ij  0,  while  no  sensor  report  was  received  when  ij«0,  l<j<k. 
The  track  history  "a"  at  scan  k  is  obtained  from  the  track  history 
"al"  at  scan  k-1  by  selecting  one  of  1+Nk  measurements  and 
incorporating  it  into  the  measurement  set  specified  by  "a".  In 
notational  terms,  a={al,  I}.  Clearly  one  history  "al"  at  scan  k-1 
gives  rise  to  1+Nk  histories  "a"  at  scan  k.  Therefore,  the  total 
number  of  hypothesis  (or  histories)  LK  at  scan  k  is  given  by 
( 1+Nk  )  *Lk_1 . 


3.2  The  Tracking  Filter  Equations 


The  equations  for  target  state  estimation  (the  tracking 
equations)  that  are  appropriate  for  use  with  MHT  data  association, 
are  listed  below.  Section  1.5  above  provides  the  notation.  Appendix 
B  explains  the  important  hypothesis  probability  equation  (11).  The 
other  equations  come  largely  from  Kalman  filtering  theory,  and  are 
very  similar  to  the  standard  tracking  equations  of  section  1.4 
above . 


Pr(k)  «  Hypothesis  probability  at  scan  k 


-  (1/0  (Pd)Nc  (i-PD)r,t-Nc  (Bft)Nf  (Bkt)Nn 

*  nj=1Nc  n[ yk-Hkxek |k_j , s ] 

*  Pr (k-1 )  ( 11 ) 
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Ka 


k  |  k-l  “  ^a^-k|k-l  +  Pa^k|k-lHk  Salk|k_1[yk  Yek|k-1  ~  ®kKa^k|k-ll  (12) 


where 


S 


-l 


k  |  k-l 


®kPk  |  k-l^k  +  Pk 


Pak|k-1  “ 

Palk|k-1  Palk(k~lHk  ^  k  |-lHkPalk  |  k-l 

(13) 

*ak  |  k-l 

“  Kalk  |k-l 

if  no  measurement  on  the 
target  was  received 

(14) 

Pak  |  k-l 

“  Pa  1  k  |  k-l 

in  the  gate. 

(15) 

K  -  E.-^lPrdOKa^.J 

-  optimal  tracking  filter  correction 

(16) 

vector  at  scan 

Kak|k  “  Kak|k-1  -  Ak 

(17) 

Pak|k  “  pak|k-l 

(18) 

state  and  covariance  update: 

Xek|k  =  Xek|k-3  +  Ak 

(19) 

Pk|k  15  ^««i  (  Pr  (  k  ){  Pak  |k_j  +Kak|k_1Kalk|k_1]  ] 

-  AkAkT  (20) 

scan  prediction: 

xek+l|k  “  ^kXek|k 

(21) 

Pk+l|k  “  *kPk|k*k  +  ^kQk^k 

(22) 

^ak+l|k  “  *k^ak|k 

(23) 

Pak  +  l|k  “  ^kPak|k^k  +  ^kQk^k 

(24) 
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3.3  MHT  Algorithm 


This  section  outlines  the  basic  algorithm  for 
implementation  of  the  Multiple  Hypothesis  Testing  method  for  multi¬ 
target  tracking,  either  in  a  simulation  environment  or  in  real 
time.  First  the  algorithm  is  described  verbally,  followed  by  a  top 
level  flowchart. 

For  notations  and  definitions,  see  section  1.5.  In  an 
actual  implementation  the  simulated  quantities  are  received  from 
sensors . 


STEP  1:  Set  scan  counter  k«l  and  all  the  control  parameters. 

STEP  2:  Initialize  the  track  and  system  parameters; 

$  r  Q,  h,  Ri  Pal1|0.  Kal^jy,  Fr(0),  xe ^  1 0 ,  P^  |  q 
STEP  3:  Simulate  target  trajectories. 

STEP  4  Simulate  measurements. 

STEP  5:  For  each  history  (or  hypothesis)  "al",  1  <  al  <  Lk_! , 

receive  unconditional  estimation  parameters  (from  STEP  11): 

Pa-'-k|k-l'  Kalk|k-1'  Pr(k-l),  Xek|k_1, 

STEP  6:  Form  gate  about 

yek|k-l  =  H  xek|k-l 

and  select  the  Nk  sensor  reports  to  be  used  in  filter 
updating.  A  measurement  is  valid  if 

(yk  -  yekik-i)T  s‘\  iyk  -  yeMk-ii  <  g2 

(this  is  a  "g-a"  elliptical  gate). 


STEP  7:  (1)  Identify  new  histories  (or  hypotheses)  "a",  for 

1  <  a  <  (l+Nk)*Lk_1. 

(2)  For  each  hypothesis,  compute  Pr(k)  using  equation  (11). 
(see  Appendix  B  for  details.) 

STEP  8:  For  each  remaining  history  (or  hypothesis)  "a",  and  each 
track  within  "a" 

Compute  Kak|k_i  using  equation  (12). 

Compute  Pa^k-i  using  equation  (13). 
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STEP  9:  For  each  history  (or  hypothesis)  "a", 

Calculate  Ak  using  equation  (16). 

Calculate  Kak)k  using  equation  (17). 

Calculate  Pak)k  using  equation  (18). 

STEP  10:  Compute  new  optimal  updated  estimates  xek|k  using 

eq.  (19).  Compute  error  covariance  Pk)k  using  equation 

(20). 

STEP  11:  Predict  for  next  scan  data: 

Compute  Kak+1(k  using  equation  (23). 

Compute  Pak+1(k  using  equation  (24). 

Compute  xek+i|k  using  equation  (21). 

Compute  Pk+1|k  using  equation  (22). 

STEP  12:  GO  TO  STEP  5  for  next  scan  data  and/or  output  the  results. 


Figure  1  below  provides  a  top  level  flow  chart  of  this  MHT 
process,  showing  the  sequence  of  events.  The  important  function  of 
"pruning"  the  hypothesis  tree  (to  reduce  the  number  of  hypothesis 
to  consider)  is  described  very  briefly  in  appendix  C.  The 
clustering  process,  also  to  reduce  the  number  of  hypothesis  to 
consider,  is  described  in  appendix  D. 
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Initialize  System  Parameters: 
X,  P,  O,  Q,  H,  k 


For  Simulation: 
Print  &  Plot  Results 


Figure  1.  Multiple  Hypothesis  Test  Flowchart. 


4.0  JOINT  PROBABILISTIC  DATA  ASSOCIATION 


4.1  introduction 


In  this  section,  a  JPDA  algorithm  for  multiple  targets  in 
Poisson  clutter  is  presented.  This  all-neighbors  approach 
incorporates  all  observations  within  a  validation  region  about  the 
predicted  track  position  into  the  update  of  that  track.  Also,  a 
given  observation  can  be  used  to  update  multiple  tracks.  Track 
update  is  then  based  on  a  probabilistically  determined  weighted  sum 
of  all  observations  within  the  validation  region.  In  fact,  the 
method  performs  an  averaging  over  observation-to-track  data 
association  hypothesis  that  have  roughly  comparable  likelihood. 
Further,  this  is  a  target-oriented  approach,  in  the  sense  that  a 
set  of  established  targets  is  used  to  form  gates  in  the  measurement 
space  and  to  compute  posterior  probabilities,  in  contrast  to  the 
measurement-oriented  algorithm  such  as  MHT,  where  each  measurement 
is  considered  in  turn  and  hypothesized  to  have  come  from  some 
established  track,  a  new  target,  or  clutter  (false  alarm). 


4  .  Z  A  J  PDA  Aly  U  L  A  tlim 


Step  1:  Set  k=l  (time  index)  and  all  the  control  parameters. 

Step  2:  Initialize  the  system  Darameters;  $,  Q,  H,  R. 

Step  3:  Either  simulate  the  target  trajectories  and  sensor 
measurements  (including  clutter  etc.),  or  in  a  real 
implementation,  receive  and  register  (scale,  time  tag  etc.) 
the  sensor  data. 


Step  4:  Validate  the  measurements  using  "g-a  ellipsoid"  gating 

test  (eq.(10)).  Form  the  validation  matrix  8(#):  (examples 
can  be  found  in  [13]). 


2jt  - 


j=l ,  2  ,  .  .  .  ,Nm;  t«=l ,  2  ,  .  .  .  ,  Nt 


(25) 


where 

w(#jt)  =  1 


0 


if  meas.  j  is  within  the  gate  of 
target  t  (event  occurs)  (26) 

otherwise 


Step  5:  For  each  t=l,2,...,Nt  compute  the  residual: 

*jt  -  Y(j't)  -  yek 


(27) 


Step  6:  Compute  covariance  of  Ylt  at  k: 
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(28) 


S(j„t)  “  Hk  Pkjk_1  Hk  +  Rj(  • 

Step  7:  Using  Bayes'  formula,  compute  the  joint  event  prob. 


Phyp(l)  -  Pr{#jzk} 

«=  ( c*/c )  *  nj=1Nt(;i’ [  -o  .5*yTjts_1yjt  /  /{ ( 2n)Nm|  s(  j  ,t)  | }  ) 

*  nt=1Nc(t,[pDt} 

*  nt„1"c,,t,U-pDt]  (29) 

where  S-1  -  inverse  of  S(j,t)  at  scan  k  (30) 

and  where  the  numbers  Nt(j),  Nc(t)  and  Nc'(t)  can  be 
calculated  from  the  validation  matrix  S2. 

Step  8:  For  each  target  t=l , 2 , . . . , Nt , 
compute 


B(j,t)  =  li  =  1Kct[Phyp(i)*w(  j,t?i)  ] 

and  for  t=0,  compute 

B  (  0  ,  t )  =  1  ~  Zj  =  1Nn,(B(  j  ,t)  ] 

Step  9:  For  each  t=l,2,...,  Nt, 
compute 

Yt  -  E3  =  1Nm(B(  j,t)*Yjt]  . 

Step  10:  Compute  the  filter  gain  AAk : 

AAk  -  HkT 


Step  11:  Update  the  state  and  covariance; 
for  each  t*l , 2 , . . . ,Nt , 


xek|k  *  *ek|k-i  +  AAk**t 

Bk  1 1;  =  Pk|k-i  “  AAk  S(j,t)  AAk 


(31) 

(  32  ) 

(  33) 

(34) 

(35) 

(36) 


Step  12:  Predict  for  the  next  scan  data: 

xek+i  ik  -  *k  xek  |  k  (37) 

Pk+i|k  =  *kPk|k^kT  +  GkQkGkT  (38) 

k+i  xe 


Yek  +  1|k  «  H 


k  +  l  |  k 
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(39) 


Step  13:  If  desired,  calculate  and  output  the  estimated  current 

target  position,  velocity  etc.  (at  time  k)  from  the  current 
state  estimate,  xekik,  or  the  predicted  target  position 
etc.,  (at  time  k+1 )  using  the  predicted  state  estimate 
xek+1|k.  The  associated  covariance  may  also  be  desired. 

Step  14:  Set  k-k+1  and  go  to  the  Step  3. 


Remark ; 


In  the  above  algorithm,  the  joint  event  probabilities  are 
computed  using  equation  (29)  which  is  derived  under  the  assumption 
that  the  probability  mass  function  (PMF)  of  the  false  measurements 
is  given  by  the  Poisson  PMF  (see  equation  (9~42)  of  [2]): 

//(m)  ■  exp{-Xvk}  (  XVk  )*  /  mi  (40) 

where  X  is  the  spatial  density  of  false  measurements ( i . e . ,  the 
average  number  per  unit  volume)  and  Vk  is  the  volume  of  the 
validation  region.  Thus  XVk  is  the  expected  number  of  false 
measurements  in  the  gate  and  m  is  the  number  of  false  measurements. 

If  the  PMF  is  uniformly  distributed,  then  the  expression 
for  the  probability  should  be  modified  accordinaly  (see  equation 
(9.46)  of  [2] )  . 
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5.0  CONCLUSIONS 


In  this  report  three  quite  different  algorithms  for 
multiple  target  tracking  in  a  clutter  environment  are  presented: 
Nearest  Neighbour  Standard  Filter  ( NNSF ) ,  Multiple  Hypothesis  Test 
(MHT)  and  Joint  Probabilistic  Data  Association  (JPDA).  For  further 
development  towards  implementation  of  the  algorithm,  the  following 
comments  and  suggestions  are  offered. 


Comments.; 


1)  The  main  dravback  of  the  NNSF  algorithm  is  that  the  tracking 
performance  or  accuracy  of  the  filter  may  become  very  poor  in  a 
dense  target  environment  because  of  possible  misassociation 
(choosing  an  incorrect  measurement  for  processing  by  a  target 
tracking  filter).  It's  main  advantage  is  that  it  is  easy  to 
implement  and  computationally  feasible. 

2)  In  the  MHT  method,  the  main  drawback  of  the  algorithm  is  that 
in  a  dense  target  environment  the  number  of  hypotheses  can  increase 
exponentially  with  each  scan,  leading  to  computational  burden 
problems.  Hence,  for  implementation,  the  development  of  an 
efficient  way  of  pruning  the  hypotheses  tree  is  necessary. 

3)  Another  disadvantage  of  the  MHT  method  is  that  the  data 
association  decision  is  often  deferred,  and  thus  a  single  best 
estimate  for  the  target  tracking  solution  is  not  always  available. 
If  target  execution  is  desired  before  the  MHT  solution  is  resolved, 
this  can  be  a  problem.  This  problem  can  likely  be  overcome  by 
providing  a  NNSF-type  solution  at  all  times  (using  the  most 
probable  current  hypothesis:  ie.  using  all  measurements  up  to  the 
current  time ) . 

4)  One  advantage  of  the  MHT  method  is  that  it  provides  a  systematic 
track  initiation  procedure. 

5)  Another  advantage  of  the  MHT  method  is  that  it  is  most  likely 
to  have  the  correct  association  solution  as  one  of  it's  hypothesis 
(hopefully  the  hypothesis  assigned  highest  probability). 

6)  The  JPDA  algorithm  is  a  non-back-scan  (or  zero-scan)  approach, 
meaning  that  all  hypotheses  are  combined  after  computation  of  the 
probabilities,  for  each  target  at  each  time  step.  One  problem  with 
this  method  is  that  it  implies  that  incorrect  measurements  are 
routinely  (and  purposely)  used  by  the  tracking  filters,  albeit 
hopefully  with  a  lower  weighting  than  the  correct  measurement.  This 
method  does  however  have  the  advantage  of  being  more 
computationally  efficient  than  an  n-scan  algorithm  in  a  heavy 
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clutter  environment  such  as  sonar  tracking. 

7)  Another  difficulty  with  the  JPDA  method  is  in  the  implementation 
of  track  initiation.  Since  the  MHT  algorithm  provides  a  systematic 
track  initiation  procedure,  an  effective  way  of  combining  these  two 
methods  may  be  of  interest. 


Suggestions : 


1)  Comparisons  of  the  accuracy  performance  of  these  different 
algorithms  for  various  target  dynamic  and  clutter  models. 

2)  Comparison  of  the  computational  burdens  of  these  algorithms  for 
various  target  densities  and  clutter  models. 

3)  investigation  of  ways  and  means  of  overcoming  the  shor tcommings 
of  these  methods,  especially  pruning  for  the  MHT  method,  and 
providing  real  time  output  from  the  MHT  method. 
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APPENDIX  A:  HYPOTHESIS  TREE  AND  HYPOTHESIS  MATRIX 


Let  2k  be  the  set.  of  association  hypotheses  up  to  time  k. 
This  set  is  obtained  from  Qk_1  and  the  latest  set  of  measurements 
Yk  as  follows.  New  hypotheses  2k-1  first  measurement  ylk,  then 
augmenting  the  resulting  set  by  associating  y2k,  etc.  The  possible 
association  for  the  i-th  measurement  yik  are 

a.  It  is  the  continuation  of  a  previous  history  (or  track). 

b.  It  is  a  new  track  (or  target). 

c.  It  is  a  false  alarm  (or  clutter). 

Each  target  can  be  associated  with  at  most  one  current 
measurement,  which  has  to  fall  in  its  validation  region. 


Example : 

For  the  configuration  of  targets  and  measurements  shown  in 
Fig.l,  the  hypothesis  tree  are  formed  using  the  rule  mentioned 
above.  In  the  hypothesis  tree  shown  in  Fig. 2,  each  node  represents 
the  track  qualities: 

"0"  is  the  false  target  or  false  alarm. 

"1,2,..."  are  the  confirmed  or  new  targets. 

The  following  figure  2  shows  the  hypothesis  tree  representation  of 
the  hypothesis  formation  technique  outlined  above-  Each  node  of 

further  branches  are 
is  considered. 


Measurements . 


the  tree  represents  an  alternative  hypothesis; 
added  to  each  node  as  a  new  measurement  point 


Measurement 


Figure  Al .  Configuration  of  Targets  and 
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after  measurement  13 
after  measurement  12 
after  measurement  11 


hypothesis  numbers 


prior 

hypothesis 


(HI) 

( Hi  2  ) 
(H18) 

(H5) 
(H22  ) 

(H8) 

(H15) 

(H25) 

<H2) 

(  Hi  3  ) 
(H19) 

(H6) 

(H23) 

(H9) 
(H16) 
(H26  ) 

(H3) 

(H20) 

(  H10  ) 
(H27) 

(H4) 

(  Hi  4  ) 
(H21) 

(  H7  ) 

(  H  2  4  ) 

(Hll) 
(  Hi 7  ) 
(H28) 


|  0:  false  target  | 
j  1,  2:  target  1  and  2  | 
j  3,4f5:new  targets  j 


Figure  A2 .  Hypothesis  Tree  for  Configuration  in  Figure  A1 . 
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The  hypothesis  matrix  corresponding  to  the  tree  can  be 
formed  as  follows. 


0 

1 

2 

3 


I 

I 

0 

0 

0 

0 


Hi 
H2 
H  3 
H4 


0 

1 

3 

0 

1 

2 

3 


H5 
K6 
H7 
H8 
H9 
Hi  0 
Hi  1 


after 

the 

measurement 

11 

after 

the 

measurement 

12 

after 

the 

measurement 

13 

002  H12 

1  0  2  H13 

3  0  2  Hi  4 

0  4  2  H15 

1  4  2  H16 

3  4  2  H17 

0  0  5  Hi  8 

10  5  Hi  9 

2  0  5  H2C 

3  0  5  H21 

0  2  5  K22 

125  H23 

3  2  5  H24 

045  H25 

1  4  5  H26 

245  H27 

345  H28 

I 

| _  hypothesis  numbers 


Figure  A3.  Hypothesis  matrix  corresponding  to  the  hypothesis  tree. 
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APPENDIX  B:  HYPOTHESIS  PROBABILITIES 


In  the  equation  for  hypothesis  probability  in  section  3.2 
above  (eq.  (11)),  an  expression  suggested  by  Reid  (eq.  (16)  of  ref. 
[11])  is  used.  Further,  the  formula  is  valid  for  type  1  sensors 
such  as  a  scanning  radar  (which  provide  a  "complete"  picture  of 
some  coverage  area  in  each  data  set,  and  hence  can  provide  number- 
of-target  information).  If  type  2  sensors  such  as  ESM  or  a  tracking 
radar  are  assumed  (which  in  each  report  only  provide  information 
about  individual  targets),  some  modifications  have  to  be  done  as 
given  in  [11].  Similar  expressions  for  computation  of  hypothesis 
probabilities  in  the  non-recursive  form  can  be  found  in  Bar-Shalom 
[2]  . 

The  equation  (11)  can  be  easily  implemented  since  it  is 
given  by  a  recursive  form.  If  all  the  prior  hypotheses  are  first 
multiplied  by  (l-PD)Nt,  then  as  a  branch  is  created  for  each 
measurement  and  its  hypothesized  origin,  the  likelihood  of  the 
branch  is  found  by  multiplying  the  prior  probability  by  either  Bft, 
Bkt,  or 

pD*N[yk~Hkxek|k_1 ,  S]/(  1-PD) 

as  appropriate.  After  all  such  branches  are  generated  the 
likelihoods  are  then  normalized. 


Example 


In  Figure  B2,  the  hypothesis  probabilities  are  computed  for 
the  hypothesis  tree  shown  in  figure  Bl,  in  terms  of  its  prior 
hypothesis  probability  and  the  probability  of  each  branch. 


s 


Pr(Hll)  -  Pr ( Hi  )*Pr(Tl) 
Pr  (  H12  )  -  Pr  (  Hl.l )  *Pr  (  T2  ) 
P r ( H 1 3 )  -  Pr ( H12 ) *Pr ( T3 ) 
Pr { H14 )  -  Pr ( H12 ) *Pr ( T4 ) 
Pr ( Hi 5 )  -  Pr ( Hll ) *Pr ( T5 ) 
Pr { H16 )  -  Pr ( H15 ) *Pr ( T6 ) 
Pr ( H17 )  «  Pr ( H15 ) *Pr { T7 ) 


Prob.  of  each  track. 

Prob.  of  prior  hypothesis. 
Prob.  after  new  measurements . 


Figure  B2 .  Hypothesis  Tree  Branch  Probabilities. 


In  figure  B2  the  probability  of  each  track  is  calculated 
using  the  simple  rule  in  accordance  with  the  track  quality: 


1/  If  the  measurement  belongs  to  false  alarm  (or  false 
target),  the  probability  of  the  track(or  br  anch )  is  given 
by 

(Bft)*(l— PD)/c. 

2/  If  the  measurement  "j"  belongs  to  i-t.h  track,  the 
corresponding  probability  is  given  by  PD*N[ 0 , S ( i ) ] ,  where 
N[a,b]  is  the  Normal  (Gaussian)  density  function  associated 
with  the  assignment  of  the  j-th  observation  to  the  i-th 
track  of  the  prior  history  and  is  defined  by 

N[  0 ,  S (  i  )  ]  -  exp(  -dj  .;*d;  </2  ) 

/{ (2n)  |S(i')  |  }  (Bl) 

for  the  residual  matrix  S(i)  and  normalized  distance 
function  di;j. 

3/  If  an  observation  belongs  to  new  target,  the  corresponding 
probability  is  computed  by 

Bkt*(l-P„)/c. 


Remark : 


If  the  track  deletion  option  is  included ,  then  the  factor 
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e~1/D*  should  be 
calculations  (as  per 


properly 

reference 


introduced  in 
[4]  pp.  255-260). 


the 


probabili ty 
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Since,  in  the  worst  case,  the  number  of  hypotheses  can  grow 
exponentially  with  time,  there  is  a  clear  need  to  limit  this 
number.  Some  schemes  for  doing  so  include: 

a.  The  first  opportunity  for  limiting  hypotheses  is 
to  require  an  observation  to  satisfy  a  gating 
relationship  before  any  of  the  possible  track 
associations  are  to  be  considered  to  be  potentially 
valid . 

b.  The  JPDA  algorithm,  by  computing  all  the 
measurements  at  the  current  time,  considers  only  the 
number  of  known  targets  with  a  single  hypothesis  per 
target . 

c.  The  "N-scan-back"  concept,  by  combining  all 
histories  that  have  common  measurements  from  k-N  to 
k . 

d.  An  alternative  method  is  to  combine  those 
hypotheses  that  similar  effects,  i.e.,  same  number 
of  targets  but  with  slightly  different  state 
estimates.  The  mean  and  covariance  of  the  resulting 
estimate  is  a  combination  of  the  individual  estimates 
and  covariance.  At  the  same  time,  hypotheses  with 
negligible  probabilities  are  eliminated. 


The  method  c  has  been  found  to  be  most  effective  in 
practice,  and  it  provides  an  efficient  approximation  of  the 
method  d.  The  process  of  pruning  hypotheses  is  highly 
dependent  on  the  applications.  Typical  ways  of  pruning  can  be 
summarized  as: 

1)  To  remove  hypotheses  with  probabilities  that  fall 
below  some  predetermined  threshold.  Disadvantage  of 
this  method  is  that  it  does  not  take  into  account  the 
computational  capacity. 

2)  To  allow  only  a  predetermined  number,  say  M,  of 
hypotheses  to  be  maintained  by  ranking  the  hypotheses 
and  choosing  only  the  M  most  likely  ones,  as  measured 
either  by  the  probabilities  or  the  score  functions. 

3)  To  rank  and  sum  the  probabilities  of  the  more 
likely  hypotheses.  When  this  sum  exceeds  a  threshold 
the  remaining  hypotheses  are  then  eliminated. 
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A  cluster  is  a  group  of  hypotheses  and  associated  tracks 
that  does  not  interact  with  any  other  group  of  hypotheses 
(contained  within  other  clusters).  The  hypotheses  within  a  cluster 
will  not  share  observations  with  the  hypotheses  within  any  other 
cluster.  The  basic  purpose  of  clustering  is  to  divide  the  large 
tracking  problem  into  a  number  of  smaller  ones  that  can  be  solved 
independently.  This  can  greatly  reduce  the  number  of  hypotheses 
that  must  be  maintained. 

The  steps  in  clustering  are  as  follows.  Initially,  one 
cluster  is  set  up  for  each  confirmed  target.  Each  new  measurement 
is  associated  with  a  cluster  if  it  falls  in  the  validation  region 
of  any  track  from  that  cluster.  A  new  cluster  is  initiated  any  time 
an  observation  is  received  that  does  not  fall  within  the  gates  of 
any  track  contained  in  an  existing  cluster.  The  cluster  is 
initiated  on  the  observation  using  the  alternatives  (true  target 
or  false  alarm)  associated  with  its  source.  In  order  that  clusters 
remain  distinct,  the  gates  of  the  tracks  within  the  clusters  must 
not  overlap.  Thus,  when  an  observation  falls  within  the  gates  of 
two  or  more  tracks  from  different  clusters,  the  clusters  are 
merged.  The  merging  must  be  done  before  the  observation  is 
processed.  If  an  observation  is  associated  with  mere  than  one 
cluster ,  then  those  clusters  are  combined  into  a  super-cluster. 
If  tracks  within  a  cluster  separate  spatially  and  have  no  more 
common  measurements,  the  corresponding  cluster  is  subdivided 
accordingly  into  smaller  clusters  that  can  be  handled  separately. 
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